基于全基因组测序的蒙古马运动相关的正向选择基因筛选

doi:10.11843/j.issn.0366-6964.2019.10.004

摘要/Abstract

摘要： 旨在从基因组水平对蒙古马运动性状涉及的相关基因进行分析。本研究利用全基因组测序技术，对蒙古马的血液进行采样和基因组测序，并且结合人、牛、小鼠等7种哺乳动物的基因组数据，对蒙古马的特有基因和正向选择基因进行筛选。分析数据得到已注释的蛋白编码基因19 475种，并筛选出了38种蒙古马特有基因和45种正向选择基因。根据GO富集、KEGG通路富集分析和前人研究，发现包括MAPK13、CCN2、FSCN3、BTG1、TF等17种正向选择的蛋白编码基因可能参与蒙古马的运动过程。这些运动相关正向选择基因的生物学功能涉及到肌动蛋白丝结构、心肌收缩的正调节、钙离子结合、调节细胞内pH、血管生成、胶原生物合成过程的正调节等。而7种miRNAs（包括eca-miR-302c、eca-miR-8974、eca-miR-9077等）则通过与mRNA的3'UTR区域作用分别潜在靶向调控运动相关的MAPK13、LOC100057755、EHD4、MICAL3、TMBIM1和F11R。以上涉及的运动相关正向选择的蛋白编码基因可以促进人们对于蒙古马运动性状的了解，并为靶向竞技型的家畜育种提供分子理论基础。

Abstract: The purpose of the work was to explore the athletics-related genes of Mongolian horse at genome level. Based on whole-genome sequencing technology, the blood of Mongolian horse was sampled and sequenced. The specific genes and positive selection genes of Mongolian horse were further detected combined with the mammalian genome databases of 7 species including human, bovine, mouse, etc.. The 19 475 annotated protein-coding genes were obtained via data analysis, while 38 specific genes and 45 positive selection genes were screened in Mongolian horse. By GO enrichment, KEGG pathway enrichment analyses and previous studies, the 17 protein-coding genes of positive selection including MAPK13, CCN2, FSCN3, BTG1 and TF, etc. were considered to be probably involved in athletic performance of Mongolian horse. The biological functions of the above 17 protein-coding genes contained actin filament organization, positive regulation of cardiac muscle contraction, calcium ion binding, regulation of intracellular pH, angiogenesis, positive regulation of collagen biosynthetic process, etc.. The mRNA 3'UTR of MAPK13, LOC100057755, EHD4, MICAL3, TMBIM1 and F11R were targeted by 7 miRNAs (such as eca-miR-302c, eca-miR-8974, eca-miR-9077, etc.). The above athletics-related protein-coding genes of positive selection can make people comprehend athletics character of Mongolian horse, while also provide molecular theoretical basis for targeting competitive livestock breeding.

中图分类号:

S821.2

潘静, 吴凯峰, 周欢敏, 张焱如. 基于全基因组测序的蒙古马运动相关的正向选择基因筛选[J]. 畜牧兽医学报, 2019, 50(10): 1985-1996.

PAN Jing, WU Kaifeng, ZHOU Huanmin, ZHANG Yanru. Athletics-related Positive Selective Genes Screened in Mongolian Horse by Whole-genome Sequencing[J]. ACTA VETERINARIA ET ZOOTECHNICA SINICA, 2019, 50(10): 1985-1996.

参考文献 34

[1]	芒来,杨虹.蒙古马遗传多样性研究进展[J].遗传,2008,30(3):269-276.
	DUGARJAVⅡN M,YANG H.Progress in the study of genetic diversity of Mongolian horse[J].Hereditas (Beijing), 2008,30(3):269-276.(in Chinese)
[2]	栾桂芝,郭聪.蒙古族三竞技与现代体育三竞技的比较及其发展[J].大连民族大学学报,2007,9(6):140-141.
	LUAN G Z,GUO C.Mongolian three-game competition and modern three-game competition[J].Journal of Dalian Nationalities University,2007,9(6):140-141.(in Chinese)
[3]	ZHANG G J,COWLED C,SHI Z L,et al.Comparative analysis of bat genomes provides insight into the evolution of flight and immunity[J].Science,2013,339(6118):456-460.
[4]	QIU Q,ZHANG G J,MA T,et al.The yak genome and adaptation to life at high altitude[J].Nat Genet,2012, 44(8):946-949.
[5]	GU J J,ORR N,PARK S D,et al.A genome scan for positive selection in thoroughbred horses[J].PLoS One,2009, 4(6):e5767.
[6]	LI R Q,FAN W,TIAN G,et al.The sequence and de novo assembly of the giant panda genome[J].Nature,2010, 463(7279):311-317.
[7]	LI R Q,ZHU H M,RUAN J,et al.De novo assembly of human genomes with massively parallel short read sequencing[J]. Genome Res,2010,20(2):265-272.
[8]	HARRIS B.LASTZ:release1.02.00,built January 12,2010[EB/OL].http://www.bx.psu.edu/miller_lab/dist/README.lastz-1.02.00/README.lastz-1.02.00a.html.
[9]	GERTZ E M,YU Y K,AGARWALA R,et al.Composition-based statistics and translated nucleotide searches:improving the TBLASTN module of BLAST[J].BMC Biol,2006,4:41.
[10]	SHE R,CHU J S C,WANG K,et al.GenBlastA:enabling BLAST to identify homologous gene sequences[J].Genome Res, 2009,19(1):143-149.
[11]	BIRNEY E,CLAMP M,DURBIN R.GeneWise and genomewise[J].Genome Res,2004,14(5):988-995.
[12]	STANKE M,KELLER O,GUNDUZ I,et al.AUGUSTUS:ab initio prediction of alternative transcripts[J].Nucleic Acids Res,2006,34(Web Server issue):W435-W439.
[13]	SAFDAR M.GenScan:bioinformatics software for structure prediction and analysis of -gene[EB/OL]. (2010-07-21). https://www.biotecharticles.com/Bioinformatics-Article/GenScan-Bioinformatics-Software-For-Structure-Prediction-and-Analysis-of-Gene-284.html.
[14]	KENT W J.BLAT-the BLAST-like alignment tool[J].Genome Res,2002,12(4):656-664.
[15]	HAAS B J,SALZBERG S L,ZHU W,et al.Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments[J].Genome Biol,2008,9(1):R7.
[16]	LI H Q,STEPHAN W.Maximum-likelihood methods for detecting recent positive selection and localizing the selected site in the genome[J].Genetics,2005,171(1):377-384.
[17]	YANG Z H,BIELAWSKI J P.Statistical methods for detecting molecular adaptation[J].Trends Ecol Evol,2000, 15(12):496-503.
[18]	YOUNG M D,WAKEFIELD M J,SMYTH G K,et al.Gene ontology analysis for RNA-seq:accounting for selection bias[J].Genome Biol,2010,11(2):R14.
[19]	MAO X Z,CAI T,OLYARCHUK J G,et al.Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary[J].Bioinformatics,2005,21(19):3787-3793.
[20]	CARLSSON L,LIND B,LAAKSONEN M S,et al.Enhanced systolic myocardial function in elite endurance athletes during combined arm-and-leg exercise[J].Eur J Appl Physiol,2011,111(6):905-913.
[21]	HIEMSTRA J A,VETETO A B,LAMBERT M D,et al.Chronic low-intensity exercise attenuates cardiomyocyte contractile dysfunction and impaired adrenergic responsiveness in aortic-banded mini-swine[J].J Appl Physiol (1985),2018,124(4):1034-1044.
[22]	WISLØFF U,LOENNECHEN J P,FALCK G,et al.Increased contractility and calcium sensitivity in cardiac myocytes isolated from endurance trained rats[J].Cardiovasc Res,2001,50(3):495-508.
[23]	HUARD J.Stem cells,blood vessels,and angiogenesis as major determinants for musculoskeletal tissue repair[J].J Orthop Res,2019,37(6):1212-1220.
[24]	LI H,LI J Y,JIANG J,et al.An osteogenesis/angiogenesis-stimulation artificial ligament for anterior cruciate ligament reconstruction[J].Acta Biomater,2017,54:399-410.
[25]	WILHELM E N,GONZÁLEZ-ALONSO J,PARRIS C,et al.Exercise intensity modulates the appearance of circulating microvesicles with proangiogenic potential upon endothelial cells[J].Am J Physiol Heart Circ Physiol, 2016, 311(5):H1297-H1310.
[26]	OTA S,UEHARA K,NOZAKI M,et al.Intramuscular transplantation of muscle-derived stem cells accelerates skeletal muscle healing after contusion injury via enhancement of angiogenesis[J].Am J Sports Med,2011, 39(9):1912-1922.
[27]	IVERSEN N,KRUSTRUP P,RASMUSSEN H N,et al.Mitochondrial biogenesis and angiogenesis in skeletal muscle of the elderly[J].Exp Gerontol,2011,46(8):670-678.
[28]	AKERSTROM T,LAUB L,VEDEL K,et al.Increased skeletal muscle capillarization enhances insulin sensitivity[J]. Am J Physiol Endocrinol Metab,2014,307(12):E1105-E1116.
[29]	BOOTH F W,RUEGSEGGER G N,TOEDEBUSCH R G,et al.Endurance exercise and the regulation of skeletal muscle metabolism[J].Prog Mol Biol Transl Sci,2015,135:129-151.
[30]	PARK J Y,WANG P Y,MATSUMOTO T,et al.p53 improves aerobic exercise capacity and augments skeletal muscle mitochondrial DNA content[J].Circ Res,2009,105(7):705-712.
[31]	SU C M,WANG I C,LIU S C,et al.Hypoxia induced mitogenic factor (HIMF) triggers angiogenesis by increasing interleukin-18 production in myoblasts[J].Sci Rep,2017,7(1):7393.
[32]	KUMAR R,NEGI P S,SINGH B,et al.Cordyceps sinensis promotes exercise endurance capacity of rats by activating skeletal muscle metabolic regulators[J].J Ethnopharmacol,2011,136(1):260-266.
[33]	CAI M X,WANG Q A,LIU Z W,et al.Effects of different types of exercise on skeletal muscle atrophy,antioxidant capacity and growth factors expression following myocardial infarction[J].Life Sci,2018,213:40-49.
[34]	ALIM M A,DOMEIJ-ARVERUD E,NILSSON G,et al.Achilles tendon rupture healing is enhanced by intermittent pneumatic compression upregulating collagen type I synthesis[J].Knee Surg Sports Traumatol Arthrosc,2018, 26(7):2021-2029.